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ABSTRACT 

For decades, Hq and go were the quest of cosmology, as they promised to characterize our "world 
model" in a mo del- independent way. Using simulated data, we show that qo cannot be both accurately 
and precisely determined using distance indicators. While Ho can be both accurately and precisely 
determined, to avoid a small bias in its direct measurements (of order —5%) we demonstrate that 
Ho/Qm (assuming flatness and w = —1) is a better choice of two parameters, even if our world model 
is not precisely AC DM. We illustrate with the analysis of the Constitution set of supernovae and 
indirectly infer qo = —0.57 ± —0.04 . Finally, we show that it may be possible to directly determine 
qo using the time dependence of redshifts, a method far less susceptible to the biases that plague 
measurements using distance indicators. 

Subject headings: cosmological parameters — methods: numerical — supernovae: general 



1. INTRODUCTION 

Allan Sandage once characterized c osmology as the 
search for two numbers (Sandage1 ll970[ ): Hq, the present 
expansion rate, and qo, the present deceleration parame- 
ter. Few would argue with the statement that cosmology 
is a much grander enterprise today and that our "world 
model" is better described by a larger set of physically 
motivated parameters (including the energy densities of 
radiation, dark matter, dark energy, and the equation- 
of-state of dark energy). In this paper we show that in 
fact, the deceleration parameter cannot be both accu- 
rately and precisely measured. While the expansion rate 
can be, avoiding a biased measurement requires a better 
choice of parameters than qo and Hq, e.g., CIm and Ho- 

Sandage introduced Hq and qo to provide a model- 
independent, kinematic description of the expansion of 
the Universe; this description begins with a Taylor series 
for the cosmic scale factor R(t), 

R{t)/Ro = 1 + H (t - t ) + ^q H 2 (t - t ) 2 + ■■■ (1) 

Where H = R /Ro and qo = -{Ro/Ro)/H%. Using 
the definition of redshift, 1 + z = R a /R, and luminosity 
distance, dh = (1 + z)r(z)Ro, and the above expansion, 
the observable luminosity distance can be expressed in a 
Taylor series in redshift 

Hod L =z + ^(l-q )z 2 + 0(z 3 ) (2) 

Note too, that no assumption about the validity of 
general relativity has been made; only that space- 
time is isotropic and homogeneous and described by 
a metric theory. We note that the next order term- 



the jerk parameter 7'n- may be added (|Visserl 120041 : 
IChiba fc Nakamural[T998l) : 



R{t)/R = l+H (t-t )--q H 2 {t-t 



) 2 + - JO H 3 (t-to) 3 +- 
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H Q d L = z+-(l-qo)z 2 --(l-q -3q ! +Jo+R- 2 v H- 2 )z 3 - 



where jo = {Ro/Ro)/Ho an< ^ ^curv is the 3-curvature 
(radius) of the Universe. 

The power of the Ho/qo (or quadratic) expansion is 
that, in principle, measurements of d^{z) can be used 
to determine the present expansion rate-arguably the 
most important number in all of cosmology-and the de- 
celeration parameter. Moreover, in the simple matter- 
only cosmology of the time, General Relativity, through 
the Friedmann equations, relates go to the physical pa- 
rameter, f^o = pm/ ' Pent (pm is the matter density and 
p crit = 3H 2 /8ttG is the critical density): 
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Further, in this model, qo is related to the spatial cur- 
vature and destiny of the Universe: if go is greater than 
1/2, f2o is greater than unity and the Universe is posi- 
tively curved and will ultimately re-collapse; conversely, 
if go is less than 1/2, Q is less than 1 and the Universe 
is negatively curved and will expand forever. The case 
of go = 1/2 is the uncurved (flat) Universe that expands 
forever at an ever slowing rate. 

This is all well and good; the question we address 
here is whether or not these two (or three), model- 
independent parameters can in fact be determined with 
accuracy and precision. The answer is simple: only Hq 
can be measured with precision and accuracy. The ex- 
planation is simple as well: at low redshifts, say z < 0.2, 
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where the Taylor expansion is most accurate, poor lever- 
age on go (and jo) and peculiar velocities severely limit 
the precision; at higher redshifts, where the effect of pe- 
culiar velocities is negligible and the leverage is greater, 
the quadratic expansion does not accurately approximate 
di,{z) (see Fig. [TJ and a bias is introduced in the mea- 
surement. For qo, you can have accuracy (restrict mea- 
surements to low redshift) or precision (use high redshift 
measurements) , but not both! (We also show that includ- 
ing the cubic term in the model-independent expansion- 
jerk jo-does change this conclusion.) 

Moreover, while one might be tempted to simply aban- 
don q$ and jo given these complications, in fact they have 
found use recently to avoid a bias of order —5% that oth- 
erwise inflicts direct measurements of Hp (due to the less 
accur ate linear expansion, H d L = z) (jRiess et all 120091 
1201 ID . To better avoid such biases, we explore other two- 
or three-parameter descriptions of dh- In particular, the 
two parameters H and fijvf (where flatness and w = — 1 
are assumed), eliminates this small bias in measuring 
H and are better motivat ed. We use the Con stitution 
compilation of supernovae (Hickcn ct al. 2009b) to illus- 
trate our findings. Finally, we speculate about measur- 
ing qo by using the very small time variation of redshifts 
(cm/sec/yr). 

Our paper is organized as follows. In Section [5] we de- 
tail our distance indicator simulations and mock super- 
nova data. We present the results of these simulations 
in Section [31 including the use of the Constitution su- 
pernova dataset to provide concrete illustrations of our 
findings. Section 0] provides some discussion and conclu- 
sions. 



luminosity scatter. 



Table 1 

Cosmological models used in our simulations. 



Model 



M 



P. 



DE 



wde(z) 



90 



ACDM 
de Sitter 
open 

Einstein-de Sitter 
w-slope 



0.28 a 


0.28 
1 



0.72 a 
1 





0.28 0.72 



-lOz for z < 0.1 
-1 for z > 0.1 



-0.58 
-1 

0.14 
0.5 

0.5 



IKomatsu eFaLl l l20TlT) . 



q(z) 0.0 




Figure 2. Deceleration parameter, q(z), in the w-slope model (see 
TableQ}. 
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Figure 1. Fractional error Ad^/d^ of the quadratic expansion for 
dL in various cosmologies. Note that Afi « 2.17Adk/d£, where fi 
is the distance modulus and Ad^ and A/i are the absolute errors 
of the quadratic expansions for and /i, respectively. 



2. SIMULATIONS 

For a diverse set of cosmological models (Table [T]) we 
generate mock distance modulus and redshift (/i, z) data 
over various redshift ranges, then study the degree to 
which qo and jo are constrained when these data are an- 
alyzed with the quadratic (Ho/qo) and cubic {Ho / qo / jo) 
expansions. We then compare the degree to which Ho 
is constrained using these expansions to the constraints 
obtained when the analysis assumes ACDM. In the lat- 
ter case, Ho and CIm are the free parameters. We also 
investigate the effect of peculiar velocities and intrinsic 



2.1. Mock Distance Indicators 

We model the distance indicators as a population of 
imperfect standard candles with an intrinsic Gaussian 
absolute magnitude scatter cr int . (While we are moti- 
vated by and will eventually use type la SNe as the 
distance indicator, our results are more general.) We 
explore two possible values of tTi nt : 0.15 mag, today's 
state of the art for SNe la after light cur v e fitting (e.g. 
Kowalski et""aT1 120081; iKessler et all l2009t Ifficken et all 
2009al IRapetti et all 120071 ); and 0.02 mag, a low level 



at which peculiar velocity scatter dominates the li un- 
certainties at z < 0.2. Note the intrinsic distance modu- 
lus scatter is simply related to the luminosity scatter as 
cr int = (dfi/dL)a L = \M{a L /L). 

Each 'measured' /i is chosen from a Gaussian dis- 
tribution of width <7i n t centered on the tru e value of 
uiz). which is given by standard equations ([Weinberg] 
2008). We then apply a modest peculiar velocity scat- 
ter to the redshift data, giving the 'measured' redshift 
-l)(z pcc +l) 



1, where each z poc is chosen from a Gaus- 



sian distribution centered at zero with c r z . De c = 0-001 
(300 km/si (e.g. iRiess et al.|[200l [20071; IRapetti et al.l 
[20071 iHicken et al.l l2009ri iLampeitl et al.l 120101) . [We 

note that there is no a priori reason why peculiar ve- 
locities should be described by an independent redshift 
scatter applied to each standard candle, and in general, 
one would expect correlations due to bulk flows which 
would prevent statistical uncertainties in the the fit pa- 
rameters from decreasing as fast as 1/^/N, where N 
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is the nu mber of data points. However, previous sur- 
veys (e.g. ILampeitl et al.ll2010f) and N-body simulations 
(|Vanderveldll2008f) have suggested that such correlation 
is not significant enough to pose a large issue.] 

In parameter estimations, this peculiar velocity scat- 
ter can be incorporated as an additional uncertainty 
in \i added in quadrature to the intrinsic scatter: 
(dfi/dz)a ZtPec — 2.\7u ZyPec / z, approximating with 
the linear expansion for d^. While o thers have assumed 
the empty universe m odel here (e.g. Kcsslcr et al .1120091 : 
ILampeitl et al.l r2010f ). essentially taking the quadratic 
expansion with go — 0, we found the difference on our 
parameter estimates to be negligible. Figure [3] shows the 
overall fi scatter 07,(2), 

a^z) = ^af nt + (2.17^ pec /zy (3) 

as well as the u uncert ainties in the Constitution data 
set (|Hicken et a l. 2009b) with the intrinsic scatter cr^t = 
0.15 mag added in quadrature. 
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Figure 3. Overall distance modulus uncertainty <r M (z) (Eq. [3} 
used in our simulations with o-; nt = 0.15 (solid) and cr; nt = 0.02 
(dashed) (including intrinsic scatter and peculiar velocity scatter), 
as well as the uncertainties in the Constitution data with cr; nt = 
0.15 added in quadrature. 



2.2. Mock Distance Indicator Surveys 

For a given cosmological model, redshift range, and 
intrinsic standard candle scatter, we simulated 10,000 
surveys, each with 500 mock objects (see Section 12. 1|) 
uniformly distributed over redshift Wc use a mini- 
mum redshift of 0.01110 and explore several values of the 
maximum redshift chosen to best illustrate the compro- 
mise between accuracy at low redshift, and precision at 
high redshift. In our quadratic expansion studies, we 
use z max = 0.1,0.3,0.5. When we explore the cubic ex- 
pansion, we use z m ax = 0.7,0.9,1.5, as this expansion 
remains accurate until these higher redshifts. Redshifts 
smaller that 0.015 are unhelpful for cosmology as pe- 
culiar velocity uncertainties are too crippling for such 

1 In theory, there is a great deal of astrophysics and cosmology at 
play in determining the redshift dependence of number density of 
a cla ss of distance indicators . In pr actice, observational challenges 
(e.g. Li, Filippcnko. & Riess 2001) often bury even those effects. 
For simplicity and to illustrate our basic results, we use a uniform 
distribution over redshift with different maximum redshifts. 

2 This is the smallest redshift present in the Constitution set. 
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nearby objects (see Fig. [3]). To eliminate simulation bi- 
ases due to mock redshifts scattering out of the redshift 
range (0.015, z max ), but not scattering into it, we actu- 
ally generate more than 500 data points over a slightly 
larger redshift range, then kept only 500 data points in 
the desired range. 

We then explore biases in precision measurements of 
Hq due to the inaccuracy of the dh expansion. To avoid 
any systematics present in high-z supernova data (such 
as evolution or dust effects), typically only nearby su- 
pernovae are used to measure Hq. For these simulations, 
we take z max = 0.1, which can be identified as the cut- 
off redshift of the low- z segment of the Constitution set 
(Fig.©. 

2.3. Parameter Estimation 

For each mock survey, we estimate model parameters 
by minimizing a \ 2 statistic. For instance, to estimate 
Hq and qo using the quadratic expansion, we minimize: 

500 / ,„ s x 2 

2/rr \ / Mi — Mexpl-HQ) <7o; Zi) \ , 

x(Ho ' qo)= t l \ — o-m) — J (4) 

Then for each cosmological model, redshift range, and 
intrinsic /x scatter, we plot 2D contours in Hq/qq param- 
eter space with 68% and 95% probability content gener- 
ated from the 10,000 parameter estimations. We treat 
the 10,000 (Hq^q) pairs as samples from a 2D Gaus- 
sian, compute the covariance matrix, then plot contours 
of that Gaussian. We have verified that the error in as- 
suming the distribution is Gaussian is not significant for 
our purposes. Each contour plot shows the probability 
distribution of (Hq,qq) pairs in parameter space given 
the intrinsic /i scatter and peculiar velocity scatter the 
data are drawn from, showing any bias in and correla- 
tion between the parameter estimates. The extents of 
the contours may be taken as error ellipses for a single 
measurement. 

To estimate go and jo using the cubic expansion, we 
minimize x 2 {Hq, qo, jo) function, then marginalize over 
Hq and proceed as discussed. 

When analyzing the Constitution data set, we pick sev- 
eral illustrative maximum redshifts, though unlike with 
our mock data sets, the number of objects is not the 
same for each redshift range. The numbers of SNe up 
to z max = {0.1,0.3,0.5,1.5} are 1141, 164, 248, 397) . 
Following standard practice (e.g . Lamp eitl et al.l 120101 : 
iKessler et all 120091 : iRapetti et al.ll2007t ). we add the in- 
trinsic standard candle scatter (0.15 mag) in quadrature 
to the published /i uncertainties. 

3. RESULTS 

Figure [4] shows i?o/<7o contour plots from analysis of 
10,000 mock ACDM simulations with the quadratic ex- 
pansion, as well as contours from analysis of Constitu- 
tion data. The trend with z max (and our basic result) 
is clearly illustrated: increasing z max increases the preci- 
sion and decreases the accuracy. Further, the precision 
at small z max is fundamentally limited by the scatter due 
to peculiar velocities. This basic trend is seen for all the 
cosmological models (except deSitter, see below). 

Shrinking the intrinsic fi scatter to 0.02 mag (the 
regime where the parameter uncertainties in the z max = 
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Figure 4. qo/Ho contours from ACDM simulations analyzed with quadratic expansion, and contours from analysis of Constitution data. 
The crosshairs indicate the true (Ho, qo) = (1, —0.58) in ACDM, where Ho is in units of its true value. The compromise between accuracy 
at low redshift, and precision at high reds hift is evident both in the simulations and in the Constitution data. Constitution Ho values are 
in units of 65 km s — 1 Mpc - 1 , as used in (|Hicken et al. 2009b). Note the expanded axes of the left-most plot. 

that the Constitution data are distributed non-uniformly 
over redshift. In fact, 80% of the data are below z — 0.8, 
and nearly all the z > 0.2 data arc not limited by the 
0.15 mag intrinsic scatter (Fig. [3]). Given these facts and 
only a single data realization (as opposed to the 10,000 
realizations the simulation plots average over), one ex- 
pects fluctuations in the z max = 1.5 (red) Constitution 
contour position roughly within the z max = 0.7 (green) 
contour. To check that this is the case, we performed 
10,000 ACDM simulations using the Constitution red- 
shifts and fi uncertainties. The resulting contour is plot- 
ted in black has go = —0.52 ± 0.14, and matches better 
the go /jo biases observed in the simulations. It appears 
slightly larger than the z max = 1.5 (red) Constitution 
contour because the actual distribution is slightly non- 
Gaussian. 

Next, we explore the small biases that creep into pre- 
cision Hq measurements due to the inaccuracy of the 
expansion. We simulate distance indicator surveys 
with z max = 0.1 (as discussed in Sec. I2.2|) in ACDM 
and in two models close to ACDM but still consistent 
with current data (ACDM may not be the exact world 
model): a model with w = —0.8 that is otherwise identi- 
cal to ACDM; and a slightly closed model with w = — 1, 
n M = 0.30, and D E = 0.75 (fi K = -0.05). We then 
estimate Hq by fitting to: 

1. Hq/VLm parameterization of ACDM (assuming flat- 
ness and w = — 1) 

2. Quadratic expansion, marginalizing over qo 

3. Quadratic expansion, fixing qo = —0.58 

4. Cubic expansion, marginalizing over go an d jo 

5. Cubic expansion, fixing qo — —0.58 and jo = 1 

The results are shown in Fig. [7J 

We note that using the linear expansion for dh yields a 
bias in H (not shown) of —4% for er int = 0.15, and —5% 
for (Ti n t = 0.02, i n good agreement with with the —3% 
bias identified bv lRiess et al.l (|2009f) using real data. Us- 
ing the quadratic or cubic expansion and marginalizing 



0.1 simulations are dominated by peculiar velocity scat- 
ter), reduces the go uncertainties by a factor of ^4, 
though the biases worsen slightly due to the altered dis- 
tribution of uncertainties as a function of z (high z data 
now weighted much more relative to low z data, see 
Fig. 13]) • We note here that while our simulations with 
Omt = 0.15 mag are intended to match state of the art 
SNe la surveys, the Constitution data set contains many 
objects with much larger uncertainties and the set is dis- 
tributed non-uniformly over redshift (Fig. [3]). Thus 
while the Constitution Hq /qo contours exhibit the same 
trends as those in the <7i n t = 0.15 simulation plot, the 
Constitution uncertainties are larger. 

Similar contour plots for other cosmological models 
(Table [lj appear in Fig. [5] Again we observe the ex- 
pected trade-off between accuracy at low z, and preci- 
sion at high z. The open and EdS models are somewhat 
better approximated by their quadratic expansions than 
is ACDM (Fig. [I}, and thus, they admit slightly more ac- 
curate go estimates. Estimates of go in the w-slope model 
are especially poor owing to its rapid evolution over red- 
shift (Fig. [5]) . The de Sitter model is worth mentioning 
because it provides a check on the reliability of the Monte 
Carlo simulations. In that model, the quadratic expan- 
sion is exact, a fact reflected in our parameter estimates 
by the absence of any biases. 

Figure |5] shows the results of including the cubic term 
when analyzing data from the ACDM model in the form 
of jo /go contour plots after marginalizing over Hq- The 
overall trend is the same as before: increasing z max in- 
creases precision and decreases accuracy. Relative to the 
quadratic expansion, the accuracy is improved and we 
thus use z max = 0.7, 0.9, 1.5, and note that over the 3 red- 
shift ranges previously considered (z max = 0.1,0,3,0.5), 
the uncertainties in the go estimates increase by a factor 
of ^5. Reducing cr int to 0.02 mag helps somewhat, but 
as before worsens the biases slightly. 

Using the entire Constitution set (z max = 1.5), we es- 
timate go = — 0.64±0.14, which deviates negatively from 
the true ACDM value of —0.58, while the simulation re- 
sult (with same z max and <r int = 0.15) of —0.45 ± 0.09 
deviates positively. To explain this, we first note again 
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Figure 5. qo/Ho contours from simulations in various cosmological models (Table [TJ analyzed with quadratic expansion. The crosshairs 
indicate the true (Ho,qo) in each model, where Ho is in units of its true value. The compromise between accuracy at low redshift, and 
precision at high redshift is visible in all models except the de Sitter model, in which the quadratic expansion is exact. Estimates of go m 
the w-slope model are especially poor owing to its rapid evolution over redshift (Fig. [2}- Note the expanded axes in the left column. 
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Figure 6. qo/jo contours from ACDM simulations analyzed with cubic expansion (after marginalizing over Ho), as well as contours from 
analysis of Constitution data. The crosshairs indicate the true (jo,9o) = (1, —0.58) in ACDM. As the accuracy of the parameter estimates 
is improved (relative to estimates using quadratic expansion), we show the results of simulations over larger redshift ranges where the 
compromise between precision and accuracy is clearer. Non-uniformly spaced redshift data and larger high-z /i uncertainties explain the 
discrepancy between the Constitution (z max = 1.5) qo estimate and the <Ti n t = 0.15 simulation result. See text for discussion. The black 
contours in the Constitution plot are the result of performing 10,000 simulations with Constitution redshifts and /i uncertainties assuming 
ACDM. Note the expanded axes of the left-most plot. 

fro m CMB anisotropy and other measurements to almost 
1% (|Komatsu et aljfeoill) . and there are aspirations to 
improve that to less than a percent (|Riess et al.l 120111 : 
ISuvu et al]l20H . H remains arguably the most impor- 
tant single number in cosmology, as it sets the age and 
size of the Universe and underpins many other measure- 
ments. 

On the other hand, not only was the minus sign in the 
definition of qo unfounded, but also qo cannot actually 
be measured! The reasons are simple: at low redshifts, 
where the dh expansion is accurate, the combination of 
peculiar velocities and the small change in d,L for differ- 
ent values of qo (lack of leverage) makes qo impossible to 
measure with any precision. At higher redshift, where 
there leverage and precision is possible, the inaccuracy 
of the Hq I qo expansion leads to bias in determining go • 
There is no 'sweet spot' at intermediate redshift that al- 
lows both accuracy and precision regardless of the world 
model. 

Moreover, the mere usage of qo as the second param- 
eter for precision measurements of Ho lead s to a non- 
negli gible bias. As we and others have shown (jRiess et al.1 
l2009f ). fitting to a linear Hubble law introduc es a system- 
atic bias of o rder —5%. One solution, used bv lRiess et al.l 
(2009. 1201 is to use the H / qo/jo expansion and by fiat 
impose the 'correct' values of qo and jo obtained with 
z ~ 1 data. However, despite the fact that Ho estimates 
using this method are no longer independent of high-z 
supernovae, this method is not robust to small changes 
in w, acquiring a ~ 1% bias if our world model actu- 
ally has w = —0.8 ( Fig. [71), within the error bars of re- 
cent measurements (Hickcn ct al. 2009b; Komats u et al.l 
1201 If ). Further, neglecting the uncertainties on q and jo 
priors underestimates the real uncertainty in Ho, provid- 
ing a false sense of accuracy. In our simulations, going 
from fixing those values to freeing them increases the Hq 
uncertainty by ~ 1%, an amount that will be important 
in future precision attempts to further constrain Ho us- 
ing this method. 
Instead and what we think is better, one can use the 



Figure 7. Hq estimates from ACDM (and models close to ACDM) 
obtained by fitting cri nt = 0.02 mock data to 5 different models 
using z max = 0.1. With <T; nt = 0.15, Hg uncertainties are a factor 
of 2-3 larger and biases are comparable. 

over the other parameters reduces these biases to order 
0.5%, and assuming values for those parameters reduces 
them further (except for the w = —0.8 model). 

Even if our world model is not ACDM, the H /Qm 
parameterization is the most robust, with biases of order 
0.01%. With cTi n t = 0.02, peculiar velocities dominate 
over \x uncertainties at these redshifts, and provide the 
limiting uncertainty in Ho- They are also the source of 
the lingering ~ 0.01% biases, as the peculiar velocity 
uncertainties are estimated using measured z, not actual 
z, so objects scattering to lower (higher) z are under- 
(over-) weighted in the likelihood analysis. 

As a final check on the reliability of our simulations, 
Fig. [5] shows the analysis of ACDM simulations and of 
Constitution data with the Hq/Qm parameterization. 
Our simulations show no biases and from this Constitu- 
tion data we determine qo = 30j\//2 — 1 = — 0.57 ± 0.04, 
with better precision and accuracy than estimates using 
the quadratic or cubic expansions. 

4. DISCUSSION 

The quest to measure Hq and qo and determine our 
world model drove cosmology for almost three decades. 
The Hubble constant has now been directly me asured to 
around five percent precision ([Riess et al.ll2011l ) , inferred 
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flow v aries with time due to the expansion of the Uni- 
verse (ISandagd Il96l ILoebl H998T) (for an overview see 
iLiske et al.ll2008l and refs. therein): 



^ = (1 + Z )H -H(z) 



-zq Ho + 0{z 2 ) 



Figure 8. Hq/Qm contours from ACDM simulations analyzed with the standard Hq/Qm parametrization, plotted with Hq/Qm contours 
from same analysis on Constitution data. The crosshairs indicate the true simulation values (Hq,Qm) = (1,0.28). Simulation Ho values 
in units of true Hq, Constitution Ho values in units of 65 km s — 1 Mpc -1 , as used in HHic kcn ct al. 2009b}. Note the absence of any large 
bias in either parameter. 

more physical two-parameter description of our cosmo- 
logical model, H /Qm- If our world model is indeed 
ACDM, then these are the only two parameters needed to 
describe cIl(z). In this case, there is essentially no bias 
in determining the Hubble constant and only peculiar 
velocities limit the precision (Fig. [7]). The same figure 
shows that even if ACDM does not exactly describe our 
world model, the biases are still very small. 

Additional parameters can be added (and may be 
needed) to completely characterize dj,: for example, w 
or wo and w a , and ftk, but given the closeness of our 
world model to ACDM, Hq/£Im is a set of two param- 
eters that affords both accuracy and precision. Fur- 
ther, we note that within ACDM additional parameters 
are needed to characterize CMB anisotropy and other 
dynamical aspects of the Universe. The standard 6- 
parameter set is fish 2 (baryon density), ft Mh 2 (total 
matter density), £Ide, ns (power-law index of density 
perturbations), r (optical depth to last scattering), and 
A 2 (over all amplitude of the spectrum of inhomogene- 
ity); (e.g. iKomatsu et aT1l2011h . and even more parame- 
ters can be added. Cosmology today is much richer than 
the two numbers that Sandage used to characterize our 
world model and the goal of cosmology. 

All this being said, there is a certain elegance and 
utility to the two parameters Hq and go, even if one 
of them is difficult to directly measure. For example, 
both can be defined independently of general relativity 
and its Friedmann equations [all that is required is the 
assump tion of isotropy and homogeneity and a metr ic 
theory (jTurner fc Riessl [20021: iShapiro fc TuTneil [2001 ] . 
Together, they characterize the most basic features of 
the expansion - rate of expansion and whether the ex- 
pansion is slowing down or speeding up. Within the pa- 
rameterizations mentioned above, go can be inferred; for 
example, with the two parameters CIm and Hq: 



where the -ffo/go expansion has been used to obtain the 
final expression. This method too is subject to bias if z 
is not small. However, the potential observable, dz/dt, 
depends upon the product of Hq and qo (meaning Hq 
must be determined separately) and is linear in redshift 
(not quadratic). Of course, this is a very challenging 
measurement-the effect is on the order of a few cm sec" 1 
decade - 1 -and a reach goal for the next generation of ex- 
tremely large optical telescopes; nonetheless we explore 
it further for directly measuring go. 

The estimate for go from this method, assuming Hq 
has been determined by other means, is given by 



est 



1 dz 
Hqz dt 



O(z) 



(5) 




i 



-0.57 ±0.04 



where the numerical value has been determined from the 
Constitution data set. 

Finally, we mention a possibility for directly measuring 
go with small bias: the time dependence of cosmological 
redshifts. The redshift of an object within the Hubble 



Figure 9. Estimate of qo (Eq. [5} obtained from dz/dt measure- 
ment. 

Fig. IH1 shows gQ St as a function of the redshift z of the 
object whose time variation is being measured, assum- 
ing that the correct cosmology is ACDM and that .Ho 
is known well. For z less than 0.1, the bias is less than 
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5% (Ago ~ +0.03), and the precision is limited only by 
that of the dz/dt measurement. (Such accuracy is pos- 
sible when using luminosity distances, but only at the 
expense of large imprecision, Ago ~ 1, cf., Fig. |U) 

Cosmology has changed dramatically since Sandage 
characterized it as the quest for two numbers. It has 
become a precision science characterized by a larger set 
of more physically motivated numbers. While, go is not 
actually measurable using luminosity distance and a hin- 
drance to accurately measuring Hq, it is nonetheless of 
interesting parameter in cosmology today-and still not 
directly measured. 
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